TECHUIOAL MEI/iOIUlTDUMS 
NATIONAL ADVISORY COMMITTEE JOH ABEONAUTIOS 



No. 749 



IMPACT BUOELIlTfl OP THIN BARS IN THE ELASTIC BANGE 
POH ANY END CONDITION 
By Josef Taul) 



L-af t f ahr t f 0 1 9 chung 
Vol. X, No. 2, July 6, 1933 
Verlag von R, Oldenljourg, Munchen und Berlin 



Washington 
July 1934 



HATIOITAL ADVISORY COUMITTBE ?0H ASROHAUTICS 



TECHNICAL MBMOHAKDUM ITO, 749 



IMPACO! SIT CKL lira 07 THIN SABS IN THE ELASTIC EA.N&S 
POS ANY END CONDITION* 
By Josef Ta-at 

SUMMAB7 



Following a qualitative discussion of the complicated 
process involved in a short-period, longitudinal force ap- 
plied to an originally not quite straight "bar, the actual 
process is sutstituted "by an idealized process for the pur- 
pose of analytical treatment. The simplifications aro: 
the assumption of an infinitely high rate of propagation 
of the elastic longitudinal waves in the har* limitation 
to slender hars* disregard of material damping and of ro- 
tatory inertia, the assumption of consistently small elas- 
tic deformations, the assumption of cross-sectional dimen- 
sions constant along the har axis* tho assumption of a 
shock-load constant in time, and the assun^tion of occon- 
tricities in ono piano. 

Then follow the mathematical principles for resolving 
the differential equation of the simplified prohlem, par- 
ticularly the developatiility of arbitrary functions with 
steady first and second and intermittently steady third 
and fourth derivatives into one convergent series, accord- 
ing to the natural functions of the homogeneous differen- 
tial equation. 

The prohlem is resolved for one type of support corre- 
sponding t.o the third Eulerian loading condition (pin-ended 
at one end, claniped at the other). The resolution for sup- 
port conditions conformahle to the first and fourth Euleri- 
an load case, may he effected in the same fashion., while 
suhstautlal simplifications are afforded for the Eoning- 
Tauh aneiilysis which treated a type of support ,cpnf ormahly. 
to the second Eulerian load case. 



It 

Stossartige Knickheanspruchung schlanker Stahe im elas- 
tischen Bereich." Luf tfahrtf orschung , July 6, 1933 i 
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The premise is a course of tbe center line of the "bar 
in unstressed condition such that precisely the fundamen- 
tal oscillation occurs during the shock, therehy ensuring 
different tar forms under different shock loads. But ap- 
propriate choice of one factor makes the mutual deviation 
of the ordinate of the "bar axis small enough so as to allow 
a mutual comparison of the tost data for different shock 
loads •vrithout appreciahlo error. Strictly speaking, the 
data are applicahle only to the chosen har shapes. As to 
the effect of a minor change in har shape, we refer to the 
appendix of the Eoning-Tauh report on the same suhject (ref- 
erence 1). 

The resolution of the differential equation -affords 
the moment formulas during the effective period of the 
shock load. After the actual shock, free oscillations oc- 
cur whose amplitudes are given through the condition (de- 
flection and rate) at the end of the actual shock (expan- 
sion in series, according to the natural functions of the 
free oscillations). An upper limit is established for the 
moments aft er the actual shock* The investigation is ef- 
fected separately for shock loads smaller and greater than 
the Eulerian Iqad for mathematical and physical reasons* 

The results are evaluated hy comparison with the mo- 
ments set up hy a static load equivalent to. the shock load. 
Tho determination of those moments for loads greater than 
the Eulerian loads involves great difficulties, however, 
with stated eccentricities, so the interpretation of the 
results ahove the Eulerian load is effected "by comparison 
of the data for the range of shock loads helow the Euler 
load» 

The results are as follows: The ratio of maximum mo- 
ment due to shock and static load is for stated eccentric- 
ities unaf feet ed.hy the magnitude of the latter* Sor short 
shock periods this ratio is less, for longer shock periods 
greater than 1, Its maximum generally occurs only after 
the actual effect of the shock load (that is, during the 
free oscillations) when the greatest deflection possihle in 
this shock period is reached du rin g the actual shock period. 
The maximum value further increases with . incr eas ed shock 
loads in two (at shock load = O) cases analyzed herein to 
3,4 (at shock loads approaching the Eulerian load). The 
excess stresses pertaining to identical shock loads due to 
the shockTlike stress decrease, on the whole with increased 
shock loads. These results refer to- shock loads, lower than 
Eulerian loaS.%. ' • - . 
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IThoroas tho moments in siiook loads below Bnlerian 
load a't'tain a maxinnim wtlch. is not exceeded during any , 
shook period, they iacrease unrestrictedly with the shock 
period at shock loads aTiove Sulerian load* 

la contrast with static loads, the Euler load may 
equally "be exceeded in shock loads when the shock period is 
adequately short* 

• ■. The results arb illustrated "by numerical examples and 
the range of validity of the approximated differential 
equation -proved thero"by. 



I. IITTEODUCTIOU 



Whereas the "behavior of structural components in pure- 
ly static stresses is comparatively yrell known, the same 
may not he said of dynamic stresises. The prohahle cause 
is that the investigation of dynamic stresse.s due to ap- 
pearance of an additional dimension, that is, the time 
factor is, a priori, more complicatod than the study of 
conditions produced undor static loads* 

As a matter of fact, stress and deformation conditions 
of structures are always due to outside interferences which 
change with respect to time. Strictly speaking, their 
study therefore always affords a dynamic prohlem. On the 
other hand, the time rate of change of these outside inter- 
feroncGs is. in many cases, so slow that the ensuing innor 
mass forces are negligibly small relative to the other 
(elastic or plastic) inside forces, with the result that 
the conditions may he accurately enough described in a sim- 
ple static investigation. This is* however, no longer per- 
missible in cases wherein the outside interactions are rap- 
idly changing. 

In bridge, and specifically in airplane design, how- 
ever, the occurrence of shock-like load effects is very 
frequent, so that the static test supplies no longer any 
reliable basis for the estimation of the strength of the 
structural parts. While high factors of safety and the 
assumption of an impact factor assures the safety of 
bridges, even without extensive dynamic studies, this con- 
dition does not prevail in airplane design, on account of 
the weight increase involved. The problem of the airplane 
designer is to ascertain the true inner and outside loads 
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and to.fi^ the dimensions to the loads as closely as pos- 
sible. 

The present paper, an attempt in this direction, 
treats the hitherto neglected prohlem, of dynamic stress, 
and specifically the behavior of an . originally not quite' 
straight 'bar in 'buckling under- shock load. The study was 
iprpiapted "by a report of ;C. Eoning On the same subject (ref- 
erence 1), and which was' translated intb German and ex- 
tended "by the writer in collaboration with the DVL and the 
;Ri jks-Studiedienst voor de Luchtvaart, Amsterdam. 

II, NOTAIIOH 



shear. 

tension. 

compression. 

Eulerian buckling load, 

support reaction perpendicular to compr,es- 
sion. 

bending monent, 

abscissa measured parallel to load P. 
abscissa measured parallel to bar axis, 
ordinate of bar axis in unstressed condition, 
"amplitude" of bar axis in unstressed condi- 

tiOJQ. 

T]x = place function. 

deflection J 

'.'H^ = time function. 

sag of point of application of compression P. 

bar length (chord length in approximations). 

chord length* 



Q kg. 

S ' ; kg, 

P = - X kg. 

Pe kg, 

A,B, kg, 



U 

X 
8 

y 

€ 



kgm, 

.m, ■ 

m, 

m, 



5 
I 
a 



m, 
m, 
m. 
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i m, distance of pulsating shear from origin of 

coordinates . 

ai snbtending slope of elastic line. 

r m, radius of curvature of elastic line. 

^ m^, , cross- sectional area. 

J . . ' 'm* , axial moment of inortia of area. 

i_ m, inert ial radius. 

E kgm"®, modulus of elasticity. 

G ■ kgra*"^, shear modulus. 

p kgs® m~* , density. 

t 8, time coordinate, 

T s> shock period. 

^ = t-T 8, tine coordinate. 

T s, period of free transverse fundamental oscil- 
lation. 

p s"^, frequency, 

9. phase shifting. 

K,X m~^, proper values. 

Z (;s, I) , core functions. 

a, see equation (Q5). 

p ratio of compression to Eulerian "buckling 

^•E ~ ' load, i.e., the reciprocal value of the 

^ static huckling strength. 

^ ratio of shock period to oscillation period 

h — m . of the free transverse fundamental oscilla- 

tion. 

c, ratio of alD solute maximum moment in the 

static and dynamic case. 
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\i, constant value (compare eqxiation (31)). 

Ajj^, nth coefficient of a series expansion. 

A,B, Cf ]3,A,S, c'l , Cs ,k]. .kg , integration constants. 
IndlceS'S 

m,nr:l,2,3. . . «. the natural figures* 



0 

E. 
I ,r 

t, 



refers to quantities occurring irith tension 
P = 0. 

quantities occurring with Eulerian Isuckling 
load, 

denote quantities to the left or right of 
i- 

place function* 
time function* 



Signs of differences: 



r 



tin 



denote the 1st, 2d, 3d, 4th derivative of 
T^jj, according to x. 



^x 

T^t» ^'t 



denote 1st, 2d derivative of ri^ , accord- 
ing to t* 



III. QUALITATIVE DISCUSSION OP PROCESS 



The qualitative description of the Tmckling stress 
due to shock load is as follows: 

An ar^bitrarily supported, originally unstressed, elas- 
tic "bar is suddenly subjected to a longitudinal force con- 
stant or varia"ble with respect to time. If tho bar axis is 
straight and the load is coincident with it, a so-called 
"dilatation wave" is produced in the bar, i.e., a longitu- 
dinal wave consisting of compressions or thinnings (or 
groups of such waves) which, as known, travel at sonic ve- 
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loclty ys/p in the material of the "bar. Arrived at the 
other end, the shock wave is at least partially reflected 
on account of the there existing discontinuities in densi- 
ty and elasticity. If the end is clamped fast, or consists 
of a .free end dipped in a medium of zero density, the re- 
flection is complete. In the further course the original 
(primary) waves superpose themselves onto the reflected 
(secfondary) waves. When the excitation of the shock waves 
is harmonic and in ahsence of material damping (i.e., zero 
viscosity of material) the superposition of tho primary and 
completely reflected secondary waves produces so-called 
standing waves, that is, waves with fixed nodes. But allow- 
ing for material damping and energy scatter at the points 
of discontinuity, there still remain apart from the stand- 
ing waves which are of lower amplitude in this case than 
"before, the so-called "advancing waves." A detailed re- 
port on this phenomenon - although limited to transverse 
waves - has heen puhlished hy Katsutade Sezawa (reference 
2). 

In unharmonlc shock excitations the type of superposi- 
tion depends on the type of excitation. The case, for ex- 
ample, of a har hit by a certain mass at a certain rate, 
and of a "bar struck hy a mass of zero velocity ("suddenly 
loaded har" ) with disregard as to viscosity and presumed 
maintenance of energy (i.e., absence of energy scatter at 
the points of discontinuity of density and elasticity) has 
heen treated ty several authors ( refer en ce3lji_ Once the 
first disturbance, moving at the rate of yE/P has cov- 
ered the length of the har, the latter is in a varlahle 
stress condition throughout its length hoth with respect 
to time and place. 

After the shock, that is, after the contact "between 
the shock mass and the har has ceased, the "bar evinces 
wave motions corresponding to the free oscillation forms. 
The initial conditions of this period of motion are con- 
tingent upon the condition existent at tho end of the ac- 
tual shock. The stress condition of the har is equally 
variahlo In time and place in the period after the shock. 

ITow, if the "bar axis is not exactly straight or if 
the shock load is eccentric, the results are "bending waves 
in addition to the longitudinal waves. The rate of motion 
of these waves is, however, contrary to that of the longi- 
tudinal waves, dependent upon the wave length. Unless the 
frequency is unduly high, it is inversely proportional to 
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the xtaive I«lig^hs; at high frequencies and when the waves 
consist ■ of ■•various groups, other laws are applicahle. With 
-vanishing wave length the rate of motion approaches that 

of . the Eayleigh waves (v = 0.9194 y'"^, reference 4), The 

rate of movement of the "bending waves is also considerahly 
lower than that of the longitudinal waves. As to reflec- 
tion and "superposition of tending waves at points of dis- 
continuities i- the same laws apply as to the longitudinal 
^Wav'SB.- ■ V ■ - 

By virtue of the fact that the longitudinal waves 
mbve faster than the tending waves, the point of origin of 
the latter is not only limited to the shocked "bar end in 
the case of the tar hit ty a longitudinal force, as is the 
case for the longitudinal wave. Bather a tending wave 
comes into "being successively at every point of the "bar, 
as'-sbon as the longitudinal wave has reached the particu- 
lar point. Each one of these tending waves spreads to 
the right and left of its point of origin to te reflected 
at rthe points of discontinuity of density and elasticity 
and to superpose themselves upon each other. But the- rate 
of movement' of these tendiiig waves is in this case unlike 
the usual TDonding waves, since thoy move in a more or less 
co'mpii'ossed part of the tar rather than in an unstressed ' 
tar. According to Sezawa (reference 2), the rate of me- 
-tioh of the tending waves decreases as the compression in 
th& tar increases. Per tuckling stress due to shock, the 
difference in rate of movement of the longitudinal and the 
transverse waves will therefore te even greater than oth- 
erwise. 

According to this qualitative discourse, the process- 
es in shoak loads in tuckling are very complicated, and to 
follow them mathematically requires various simplifying 
assumptions. 

IV. SIIvIPLIFYIira ASSUtlPTIOHS 



la The first assumption is infinitely great rate of 
movement of the longitudinal waves, or in other words, 
that the tar reaches longitudinal stress condition within 
infinitely short time. Then the tending waves emanating 
from the different points of the tar are simultaneously 
excited rather than succe_s sively , and instead of" a wave 
motion there are standing oscillations to tegin w.ith. 



r 
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This assniaption is the more applicable as the frequency of 
the shock motion is lower and the shock force greater, he- 
cause hoth factors tend to lower the rate of movement of 
the tending waves "but not that of the lon^;itudinal waves. 
Moreover, it is anticipated that the conditions of the os- 
cillating "bar corresponding to this assumption will more 
closely approach actual conditions as the time interval 
"between entry of disturhance, i.e., start of shock and end 
of shock is greater. 

On these premises the analysis is restricted to the 
effect of the bending oscillations. The stress due to the 
longitudinal force is readily obtained by dividing the mo- 
mentary shock load by the cross- sectional area. 

2, The study^jjis confined to the investigation of buck- 
ling due to shock /slonder bars. From the cited qualita- 
tive analysis it follows that the energy of the shock mass 
goes partly into bending and partly into longitudinal os- 
cillation energy of the bar. The cozaponent of the bending 
oscillation increases with the slenderness of the bar as 
seen hereinafter. 

A bar pin- ended at both sides (fig. la) is subjected 
to static compression P. Presumedly, the Taar was not 
quite straight at the start and the deviation of its axis 
from a straight line corresponds to a half sine wave. It 
can be proved that the elastic line is sinusoidal in this 
case also. With the arc length as independent variable 
(reference 5), wo have: 

y=€Bin^s (l) 



■max 

Moreover, according- to figure lb: 



= ri^ax sin ? s (2) 



ds 

The chord length follows from 



^^^■^^Usina (3) 



a = /• dx = / cos a ds (4) 



o 



In second approach cos a may be replaced by (reference 
5): 



* 
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cos a-= 1 - |i = 1 - ain|_a 



0^1 "because of (3), "by 



cos 



a 



This value of cos a written in (4) gives with due 
Allowance for (1) and (2): 



= / [1 - I (^'. + 3 e r\^^^ + n^^,) cos= E s] ds = 



a 

o 



By the same argument the chord length prior to loading 

is: 

o 4 1 

Consequently, the sag of the applied moment of the load 
is: 

■ g 

6 = ao - a = — (riliax + 2 c Timax) (5) 

The deflection in har center in this case is conform- 
able to the Koning-Tauh report (reference 1); 

-S - 1 

p 

with Pji = '^-r^--^ the Eulerian ouckling load, 
With a = (5) and (6) give: 

4 1 ( 1 - a) 

This formula yields versus I, wherehy it is 

to "be rememtered that I implicitly occurs also in Pjj. 
The curves ohtained for stated P, B, and J values are on 
the order of those shown in figure 2, As the slendernoss 
of the bar increases the displacement 5 of the applied 
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moment Increases nnich more rapidly than linear even when 
c is assumed as 'being constant. On top of that e, that 
-is, the prohahle maximum deviation of the har from straight, 
is Infrequently assumed so much greater as the "bar is thinner. 
(Muller-Breslau, for instance, figures with an initial ec- 
centricity of l/20a) 

Thus* while the sag of the point of applied moment in- 
creases with the length of the har at a much faster rate 
than linear, that due to the longitudinal compression in- 
crbases only linearly with the "bar length (Hooke's law). 
Thoroforo, the more slender the har, the greater the do- 
flection compared to the displacement of the applied mo- 
ment of the force due to compression. Assuming for sim- 
plicity that the force P ., is constant in, time, the energy 
required to deflect or compress the "bar is simply propor- 
tional to the deflection or compression or« in other 
words, the ratio of "both is actually so much higher as the 
har is more slender. 

Hence the omission of the longitudinal oscillations 
appear still more justified since the investigation is 
confined to slender hars. 

3.. The damping effect due to material viscosity is also 
ignored as its effect "will "bo rather negligihlo in the 
short period processes involved here, although it.s inclu- 
sion should cause no insuperahle difficulties. 

4. All investigations retain their validity only in 
stresses helow tho elastic limit. 

5. The deviations from straight line as well as the 
.deflections are assumed to he small quantities compared to 

tiie length of tho "bar, so that tho curvature of the elas- 
tic line may he accurately enough defined with d^r|/dx^. 
The scope of validity of this assumption is analyzed in a 
suhsequent section. 

6. We assume a har with sectional area constant along 
the har axis, with constant inertia moment, and elasticity 
modulus. 

7. The transverse (tending) oscillations are treated 
as plane prohlem, i.e., all deviations from straight line 
are presumed to lie in one plane. 



r-ji- -e. As to- tile tinio^ratB of shock, it is siniply assixmed 
tlia't tlj.0 aliock' load is a constant independent of time. 
. :':rAj-t]jo-tigli,-t]ais ■ inay not exactly correspond to the actual 
shock of landing ef an airplane or other similar shock 
leads, it is a good approach in the majority of cases* 

9* The effect of the rotatory inertia Is disregarded, 

"i.e., the term p J — deleted in the differ-' 

3 t 5 X . . 

eritial equation. 

■ V. DIFPEEENTIAL^ EQUATIOHS ASH BASIS OP HESOLUTION 

1. General it ies 



The approximate differential equastion is 

E J - X + p p -^g^ _ X C8) 

(X must he written with positive sign as tension.) The 
resolution of this equation, i.e.., the determination of the 
transverse oscillations of a har loaded under constant ten- 
sion X, whose center line in.. unstressed condition is. 
y = f (x) , is readily tractalsle for the case of a "bar 
hinged at "both ends, "by expanding y = f(x) in a series 
whose individual terms correspond to' thb natural func- 
tions of the free har oscillation. 

For the generalized case (i.'e.,'any support condition) 
the expansion of f(x) conformahly to the natural func- 
tions; of the free "bar oscillation is ahortive. It rather 
necessitates . the development of the right-hand side of (a) 
according .to ..the natural functions of the homogeneous equa- 
tion 

^ ^ $ - ^ I?' = .° ■ .- ' -^'^ 

the differential equation of the straight har stressed in 
tension X.. It possesses mutually orthogonal natural 
functions,: from which arbitrary funct.iiOins may he developed 
in convergent series. 

The exact proof of this dovelopahility follows from a 
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Hll'bert theorem (reference 6), conformable to vhlch any 
ftinction ( re.f e,r.enc0 7), according to the natural fxmctions 
of a continuous symmetrical core may "be oscpanded in a con- 
vergent series. 

fhe source-like representation of the elastic line is 
obtained with the so-called "influence function" (&reen 
function). This is the equation of the elastic line with 
X as independent variable due to a force 1 acting at point 
x= |, According to whether x~ S or x= the re- 
sults are two different terms, which are denoted by 
£t(x, ^ and £r(x,0, respectively. Maxwell's law of 
mutual displacements states that a force 1 acting at point 
I . causes the same • deflect ion at x as a force 1 acting 
at' X in i. Therefore, the influence fxmctions Kx.(x,|) 
valid for^ x^ | and Ki.(x,|) valid for i are sym- 

metrical - (the appendix gives as examples the symmetrical 
influence functions for the support conditions correspond- 
ing to the second and third Euler case) - that is. 

Kl(x, i) = KrCl.x) . 

But this influence function is exactly the core function 
of the integral equation, thus proving the developability 
according to natural functions. 

Maxwell's law adduced as proof of the validity to a 
bar compressed in constant tension and bent in arbitrary 
shear follows from the fact that the premise of Maxwell's 
law: the law of linear superposition is equally valid, 
provided the tension remains unaltered. The validity of 
the law of superposition follows from the linearity of the 
differential equation (9) , which is briefly discussed in the 
next two sections with particular reference to the corre- 
sponding natural functions and period equations for the 
support conditions characterizing the four Bulorian cases. 

2. Transverse Oscillations of a Straight Bar 

.Loaded in Tension 

a) Solution of _..the h^^ differenti al equation . - 

The particular solutions of synchronous character of (9) 
in conformity with the first assumption, are: 

= Tit 'Hx 

with T|t = f(t) and t)x = f(x). Putting 



T) = sin p t e 
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ln^(9)., yields the tji-pical equation: 

. / E J K* - X - p P = 0 (10) 

of the four roots. 




K =± /---± / -.-IL^+ Blps (ii) 

/ 2EJ- V 4 EJ 

of. (IQ): Kx =* ± X are real for real values of p,' and 

two,, .Ks X.i are, imaginary, and this holds for posi- 
-fciye X (tension) as well as negative X (compression). 
Then. the solution reads: 

1 Tlx = Ach\x+BshXx+C cos + D sin. X x (12) 

Per imaginary values of p, we obtain four real 
roots, = ± Xi and Kg = ± Xi with positive X when 

assuming that 



Xf __ ^ 

4 



EJ ^ 



(13) 



For tension loads (12) then hecomes 

Tlx = A ch Xx + B sh Xx + C chll-x + D sh Xx (14) 

and for compressive loads, 

r\:^ = A. cos Xx + B sin Xx + C cos Xx + D sin Xx (15) 

Imaginary values of p occur according to VI, 3, h 
in compression exceeding the Eulerian load, The same sec- 
tion also discloses that (13) is met so long as the com- 
pression is not in excess of twice the Eulerian load. 
Compressive stresses "beyond this double Eulerian load are 
without the scope of this report, • 

It also follows from (11) that 

+ = 7? = (16) 

whorehy (9) "becomes 

[(s J x^ - X) x^ + p p Ti^] -n^ = 0 (17) 
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From (16) follows "X = \ trhen X = 0. Without ten- 
sion .the term for t^^: thus 1)6008168 

Tlx = A- ch + B sh \x + 0 cos \x + D sin \x (18) 

and (17) reduces to 

(B J X -nt + P rit) Tlx = 0 (13) 

h) Determination of co nstant s .and derivation of th e 
period equation for various s upp ort conditions.- In the 
following the constants In (12) together with the natural 
values of \ and K . are analysed for the support condi- 
tions of the four Suler cases. 

la Euler case: one end fixed, the other free. 

Boundary conditions: r\x =0 and 'n3c'=0 for x=0 

Ti3."=0 for X =1. 

The fourth Is obtained for x = I as follows: 
I'or X = I 

Tijj* = tan a sin a r- a, 

hence the shear. 

Q = - X sin a = - X Tlx' • 
On the other hand, 

Q = - B J Tlx"* 

thus 

•Hx'" - ;|rr,x' = 0 (20) 

when X = I . 

From the "boundary conditions for x 0 follows; 

A + C = 0, consequently C = - A 

— _A. 
X B + KD = 0, " D = -^B. 
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Therefope," 

Tia: = A(ch - cosHx) + B( sh \ x - ^ ain Xx) (21) 

.A 

With X = I, the Isoundar^r conditions give: 

'nx"=-A.(\® ch cos \l)+B(?^ Bli>^l sin X l)=0 (22) 

aiixi _ _ _ _ 

. . . A(\ all sin \1>+B(X ch X-l+XX cos \ I) - 

« ^X^ XAXX sh Xl+X sin \l)] +B( \ chXl-.\ cos X l)=0 

• EO;. ;. . 

ori with allowance for (16) 

A(\X^ sh sin \l)+B(\\® ch cos Xl)=0 ' (23) 

So that A and B do not disappear, the denominator 
determinant of (32) and (23) must he zero: 

\^ ch \l+ \^ cos\'i- \® sh + sin XI .. . . 

\}? sh Xl-X^f sin XI X X^ ch Xl + X^ cos XI 

The evaluation of this determinant gives the tran- 
scendent period equation: 

2 X^X^+(?^ +t )ch XI cos XZ+XX(?f -X^ ) sh XI. sin.Xl=0 (2.4) 

which, together with (16) gives the natural values of X 
and X for the different longitudinal forces. 

Por X = 0, X = X and (24) he comes 

chXl cosXl +1 = 0. 

2. Euler case: hoth ends hinged. 

Boundary conditions: T|x.=0 and Tix"=0 for 3c=0 

rix=0 and •nx"=0 for 3?= I 

Por X =s 0 , we have: 

A + 0 = 0 . 

I hence, A = .0 = 
X^ 'A - 1? C = 0 ^ ■ ■ ■ 
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for X = I, we have: 

Tlx = B \l + H Bin'Kl = 0 (25) 

-Hx" = B Bh Xl-'X^ D Biii"\l= 0 (26) 



The condltioB. according to which the denominator de- 
terminant of these equations must disappear, gives the pe- 
riod equation: 

sh XI sin \l = 0. 
This equation is fulfilled for sin "XX = 0, that is, 

for 

A, _ _., .... 14_. 

From (25) and (26) further follows 3=0, whence 
the validity of 

Tlx = D Bin "X X (27) 

Since this equation contains only \, the proper functions 
in the second Euler case are unaffected "by the longitudinal 
force. Solution (27) is equally valid when no tension ex- 
ists. For this reason it was possihle in the seeoAd Euler 
case treated in the Koning-Tauh report (reference l) to ex- 
pand the right-hand side of (8) conformahly to the natural 
functions of the free transverse oscillations. 

S. Euler case: one end clamped, the other pin-ended. 

Boundary conditions: r]j^=0 and Tix'=0 for K=0 

•nx=0 and 'nx"=0 for 70=1 

' The conditions for x = 0 are the same as for the 
first Euler case, hence (21) is valid. 

Por X = I, wo have: 

'nx=l = A(ch Xl-cos \l )+B( sh \l- ^ sialvl )=0 (28) 
•n" 7 = A(?f ch Xl+\^ cos Xl)+S(>f sh \l+X\ sin"\l)=0 (29) 
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The cbn'dit ibn of vaxxiehing denominator determinant 
affords: 

(ch Xl-cos "Xl) (X^ sh Xl+XX sin - 

- (sh \l - = sin Xl) (X® ch Xl+V cos XI) = 0 
X 

and, after appropriate conversions, the transcendent peri- 
od equation: 

^ tan XI - th XI = 0 (30) 
X _ 

The roots of (30) give the natural values X and X for 

the different X values. 

With M. = r . (28) yields: 
A 

p, = - _cli-. .>^i_-_.£osAJL (31) 

shXl-j=?sinXl 

X 

which, in conjunction with (28) and (29) gives: 

ch XI + sh XI = 0 (32) 

cos "XI + |JL ^ sin XI = 0 (33) 

X 

For vanishing longitudinal force (X = O) (30) hecomes: 

tan XI - th XI = 0 (34) 

since X = X, 

Since this paper treats the third Euler case as an 
example of the generalized method, the determination of 
the natural values requires a more detailed discussiont 

The introduction of 

a = (35) 

E J X 

changes (16) to _ 

X^ - X^ (1 - a) (36) 

and (30) to 
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tan \l til {\l ^fl " a) = 0 (37) 

^ • • • ■;■ ■ v^~^ . 

the" roots "\ 1 of this eqtiation are ohtained as inter- 
sections of the curves 

f (X) = tan "X I 

and 



f (a.X) = —3: th (\l y 1 - a) (38) 

y 1 - a 



For a = 1, 1 (a,\) = ^. consequently. 



Cth (\l Jl - a)] 

'T" 

da 



f (a, >.),,=,, = ^al- ^— r™— = 



is a straight line. 

For a > 1 (38) "becomes: 



f (a.X) = t§.?L_L>iJ N^,-, AI 

/a - 1 



Figure 3 shows the functions f OO and f(a,X) plotted 
against XI for several a values (parameter). 

For a = 0 , 

; , 5 9 4n + 1 

Aol=— IT, tTTj •»•! T""" "T 1 « • • • 

4 4 4 

is accurate enough. 

All these values are at least within 0,1 percent accu- 
racy of 

|JLo = - 1 

8 0 that, 

ry^Q =, A(:ch -Aq x - cos^ Xq- s - sh Xq x + sin X^ x) (39) 
is also sufficiently exact. 
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In figure 4 the XI values for the fundamental and 
the first higher oscillat_ion have "been plotted against a. 
Up to ahout a = 1 the XI values show a slow increase 
with a; "beginning at a > 1 the mean rise is steepor 
and the curves coil about a mean curve (dashes) which, eX" 
pressed in formula, is 

\l (40) 

1 - V a - 1 

where n = 1 = fundamental oscillation, and n = 2 = first 
higher oscillation. Equat ion_(40) is exact for whole mul- 
tiples of tt/2. Por a = 2, Xl = «. 

4. Euler case:, hoth ends clampod. 

Boundary conditions; r\y^O and 'nj'=0 for x=0 

r|3j=0 and t^x'"^ • 

The conditions for x = 0 are, as in. the first Eulor 
caso, so that (21) is applicehlo. 

Per X = I, wo have: 

r\x=i = A(ch \l-cos X«)+B (sh \l - = sin XI) = 0 

X 

n' =A( X sh Xl+X sin'Xl)+3(X ch Xl*.X cos "^-V) = 0 .- : 
x=l 

The disappearance of the denominator determinant stip- 
ulates: 

X(ch Xl-cos ■Xl)^~(sh XI - ~ sin'Xl)(X sh Xl+X sin XI) ^ 0. 



Thus the period equation is: 

2(1 - ch XI cos XI) + ------- sh XI sin XI = 0 (41) 

X X 

The X and X values are ohtained as roots of (16) and 
(41) . 

Per X = 0, X = X, so that the period equation he- 
comesi _ 

ch XI cos XI = 1. 
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V.I.. 5PLUTI03T OF THE UNHOMOGEKBOUS DIPFERElToiiAL isqUATIOU 
■ 1.. Greneralities 

With, the "proper, functions derived in V,.2, we proceed 
in the general case hy developing for given original "bar 
form .y(x) the right side of (8) according to tho proper 
Ifianctions of (9) and then poae a corresponding oq.uatlon ' 
for the place function r\y^' in tho solution Ti = "Ht 
conformable to the natural functions of (9)« . 

In view of the comparatively complex form of th€» nat- 
ural functions together with the tedious series develop- 

ment of X — & , this method would, however, "become quite 

troublesome. AddocL to this is the fact that different 
functions would havo to he sot up for the original "bar 
form y, even if intending to allow for the prohahle cases 
only. It must therefore he understood that the investiga- 
tion - exactly like similar ststic studies - is feasible 
only for stated bar forms, and that tho conditions antici- 
pated for other forms must be arrived at by deduction from 
the obtained results. 

Under these circumstances the thought lies close to 
discuss the general method with such a bar form as to lead 
to mathematically convenient formulas. The most elementa- 
ry case is a bar form developing only the fundamental os- 
cillation during the shock period. But, excepting the sec- 
ond Euler case, the natural values and through them the 
form of the fundamental oscillation differs for different 
X (see V,2,b - figs. 3 and 4). Accordingly, different 
bar forms are necessary for different X, to insure the 
appearance of the fundamental oscillation only. On the 
other hand, an analysis which presumes variable original 
bar forms as X increases, would be unsatisfactory as it 
would, to be sure, afford a comparison of the bar under a 
stated static load with an equivalent dynamic load, while, 
however, affording no possible means of coinparing the proc- 
esses occurring under different longitudinal forces* 

However, we shall use this method for the third Euler 
case-,* which is being analyzed to explain the general moth- 



*!rhe period equation of the third Euler: caso is much more 
simple than than that of the first and- fourth Euler caso 
(see (24), (30), and (41),) thus making the mathematical 
treatment of the third case somewhat more convenient than 
tho first and third case* 
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od.V It is possible to show that with appropriate handling 
of a constant the "bar forms at which precisely the funda- 
mental oscillation occurs are so little imlike under dif- 
ferent longitudinal forces as to make comparison of the 
results obtained for different longitudinal forces possible 
jtr.idb-lajout great errors* 

As concerns the effect of slight modifications in the 
"bar form, we refer to the appendix of the report cited in 
reference 1* 

2. Choice of Interference Punction (Bar Shape) 

To insure that the "curved" bar executes the fundamen- 
tal oscillation of the straight bar the second derivative 
of its original form must be chosen conformably to the 
first natural function of the straight bar. 

Accordingly, we choose for -^-s 

B^y — 2 — ^1 ~ 

•r-^- = - CiX-x Cell ^1 x-cos \i x+|j, ( sh X,i X - =^ sin K^x ) ] 

(42) 

which up to constant Ci, agrees with the first natural 
function of (9), (Compare (21) and (31).) 

Then (8) resolves to 

K ^ 

'n='nt 'Hx^'Ht Cell x-cos \i x+\x ( sh Xi X - :=r*- Bia \iX )] 

^1 (43) 
Writing (42) and (43) in (8) and omitting the index 1, 
we have: 

[(EJ \^-X)?\f Tit+pTT^^+X c'X®] (ch X 3H-IJ. sh \ x) - 

- C(EJ "X^+X)\^ Tlt+P^ *^t+^ ^ X^](cos X x+pi= sin"X x)=0. 

X 

or., because of (16) , 

C (E J \^ + x)"x^ + p p + X e ] Tijj = 0 (44) 

!Tow we try to ascertain to which original bar form-, 
the chosen function (42) corresponds. i..,.e. , find the trend 
of the eccentricities y to ensure that the shock load X 
exactly produces the fundamental oscillation* 



lil.A.C.A. lechnlcal LEomorandum "So, 749 



23 



-Twice integrating (42) gives: 
F= - c I cli \ X + cob'X X + p, f ^ sh \ + 



+ -~ sin "X X 



) 



+ Ci X + C2 (45) 

X 

Now we can prescri'be certain conditions for defining 
Ci and Cg,. The nearest One would "be to require that the 
■boundary conditions for function y agree with those of 
the elastic lino of the stressed har. 

The conditions then would 'be: 

1, For X = 0 y = 0 

2. » X = 0 7- = 0 

dx 

a, " X = I y = 0 

dV 

4. " X = I -4=0 

With only two constants available, however, only two of 
these conditions can be met. Tnere are (4/2) = 6 possi- 
hilities. 

We select that for which conditions 1 and 2 are ful- 
filled. It is simpler in many respects than other forms. 
Subsequently, it is shown that for it the maximum moment 
nearly always occurs at the point of fixity. Added to that 
is the ready comparison to har forma under different longi- 
tudinal loads cited in VI.l. Jor this har form, it is 



and 



consequently, 

•2 



y = - e 



S-ch\x + cosXx+[A(^sh?\.x+.— sin^xj- 
"•'(^ + l) (1 + M. X x)J (46) 



34 



IT.A.C.A* Tochnlcal Memo raxidum ITo 749 



Tiio evaluation of |(46)- gives the curves shown in fig- 
ure 5. The values - ^ for x = 0,5 I, x = 0.75 I and 

X = I proportional to the ordinatos of tho har form aro 
plotted versus a. The ordinates for x = 0.25 I are too 
small to be representahle at the scale used in figure 5a 

For a = 1, which corresponds to the Euler load as 
tension, as shown elsewhere, X = 0 because of (36), and 

{J, = ^ because of (31). Thus (46) may be expressed as 



y = - € 



■^J (ch \ x-l) + (co3 X x~l)+ ^ iJL(sh X x-X x) + 



+ ^ \i. { sin \ X- \ x) 



after which we determine separately the unknown terms, 



(Ch X . - 



( ch X z 
resulting for 



_ T« sh X 
- X — 2- 



= 1. Purthormor o, according to (33) 
IJ. X = -"X cot a I 



and consequently, 



[x^Jx==o I >f Jx=< 



CO ^ 



The third term in the expression for y and through 
it value y itself becomes co . Thus the curves (fig. 5) 
have a point of discontinuity when a = 1 or, in other 
words, the original bar form is undefined for a = 1, so 
that the case a = 1 must bo excluded from the invest iga- 
t ions. 

Figure 6 gives the original bar form (corresponding 
to (45)) for a = 0 at exaggerated scale. We use it as 
original form, We now reduce the bar forms resulting from 
(46) for different a in such a manner that the maximum 
ordinates (= ordinatos for x = I) agree for all a val- 
ues with the maximum ordinates for a = 0, By virtue of 
(32) and (33), (46) gives the maximum ordinate at 



ymax = e + 1^ (1 + {J. ^ I) 
V.X ' 



(47) 



N.A.C.A. Tecliulcal Memorandum No. 749 25 



Let € = €q for a = 0; then, "because of 
it gives 

Cymax]a=o = 2 €o (l - | n) = - 5.86 (48) 

To assure the same value of ymax other values of 

a also, we put: 

5,86 £n 5,86 1 - a 



( 



13 \ l + p,\l2-a^° ^^^^ 

-^5- + lj(l + p. \ I) 

€ 

figure 7 shows the reduction factor -r— against a. 

o 

With this choice of €, (46) manifests bar forms for dif- 
ferent a which practically do not deviate from each oth- 
er, as seen in figure 8, where - v/^o x = I, x = 

^ I and X = ^ I proportional to the ordinates of the re- 
duced har form are plotted against a. The discrepancies 

of the curves for x = ^ ^ and x = ^ I from a line par- 

4 2 

allel with the axis of aTsscissas are very minute. The 
greatest percentage of discrepancy of ibhe ordinates among 
each other oliviously occurs wher. ^ < |-, which is not 

shown in the graph. In the ran£:Q of 0 <a < 1.09 , (higher 
a values should he of no practical interest) the maximum 

y 

discrepancy among the ^— ordinates amounts* for instance* 
1 

for X = ;j I, to ahout 15 percent. Admittedly, the cri- 
terion of the approach is not this figure hut the maximum 
error referred to the maximum ordinate. This, however, 

does not occur at low "but at ^ = j and in the range 

of 0 < a < 1.09 is loss than 2,5 percent, and in the 
range of 0 < a < 1.05 even less than 1.6 percent. 

This proves that har forms chosen according to (46) 
and containing the c value given in (39) show such lit- 
tle discrepancy as to render comparison of the computed 
results possible without committing any great errors. 
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S» Solution 
With a view to (35), (44) gives for -nt: 

p P rit + E J X* (1 - a) Tit = a E J (50) 
There are three distinct cases: 

(a) a < 1; (h) a > 1} (c) a « 1. 

He confine ourselves to cases (a) and (h), as the har form 
y is not defined for a » 1. (See VI, 2.) 

The physical interpretation of these _cases is as fol- 
lows: Let Pjj = Eulerian huckling load, = correspond- 
ing X value, and aji = ratio of longitudinal load to Eu- 
lerian load. Then 

* E J A-s 

or, hecause of (28) 

as = a (52) 

■ A. jj 

It is shown in TII.l that ,\jt agreos with the value 
of ~k for a = 1; (52) is ^sraphed in figure 9. For ag = 
1, it is a = 1, so that our two cases (a) and (h) may 
also he expressed as 

ajj < 1 and ajj > 1. 

Physically, this means: 

The shocic load is lower than the Eulerian load; 

The shock load is higher than the Eulerian load. 

a). Sh ock l oad lo-^e r than, E ule rian loa d.- With the as- 
sumed initial conditions, 

T|t = 0 and T^t = 0 t = 0 

the resolution of (£EC^) gives, analogous to the second Euler 
case in the Konxng-Tauh report (reference 1): 



N.A.C.A. Ischnlcal Memorandum No j. 749 



27 



T|t = i*"-"a ^ ^ *^ ^^^^ 

where 

p = yn /"^"-'^ ^ ^^^^ 

is the frequency. With index 0 denoting the quantities 
referred.to the free fundamental oscillation (a = O), we 
have: 

P = Po r~2 -/I - a (^^^ 

Ao ■ 

figure 10 shows the ratio versus aw for the fun- 

Po . 

damental oscillation. As in the second Buler case, the 
frequencies drop as the load increases and disappear alto- 
gether when a = 1. 

Thus with (43) and (53) the resolution of (8): 

T) = Ti^. Ti- = --^ — € (1 - cos p t) 
^ ^ 1 - a 

[ch \ x~cos X z+ij,(sh Xx--i^sinXx)] (56) 
At tho end of the shock t = *r, it is: 

^ = __a — e (1 - COB pT ) 
1 ~ a 

[ch "K x-cos ^ x+|jL(sh X x - ^ sin K x)] (57) 

K 

[ch \ x-cos \ x+p,(sh \ X - ^ sin X, x)] (58) 

A, . 

Tho free oscillation of the har after the shock is 
with duo allowance for (39): 

r\= 2 An sin(poii t+cp^) ( ch Xon x-cos Xq^q, _x-' 



n 



— sh Xon xl-sin Xqu x) (59)' 



The choice of time coordinates T after Jihe shock is 
such that t = 0 at tho ond of tho actual shock period, 
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that is, the start of the free oscillation. Then t = t - T 
and the start of the free oscillations are according to 
C59): 

'ti •= S sin 9ji(ch Xon 3= *• cos Xqb ^ "" 

•- sh ^on x + sin ^on x) (60) 

1^ = I -^n Pon cos 9a(°li ^on ^ °os \,n =^ " 

- sh ?\on ^ + sin "S^^ x) (61) 

These values of the deflection and rateof deflection 

at the start of the free oscillation must agree with the 

corresponding values for the termination of the actual 

shock (57) and (58). To define the constants and q)j^, 

in (59), we then have: 

S An sin 9a Tixon = l~^i:^^ " ^ 'Hx (62) 

g ^ Pon cos (Pn e p sin p T t^j^ (63) 

Herein: 

'nxon=ch \on 3C-C0S \on x-jjLonCsh Xoq, x-sin Xon ==) = 
= I [o^o^^ (1 + Hon) + e"^°^^ (1 - p-on) - 



i\on3c , » -i^on* 



- e"'^°»* (1 - Hon i) " e""'^oi^* (1 + Hon d] (64) 

is the nth natural function of the free "bar oscillation,* 
and again hy omitting index 1: 

T^j. = ch X,,Jc-cos X x+H (sh A x - :r sin X x) = 



1 r ^3c. . -Xx. . iXx,, Xi\ 

I [^e (1 + ij,)+ e (1 - h)- e . (1 - p, " 

-iXx / Xi\ 



- « 



(65) 



♦Simplification of Hon ~ *" 1 , is purposely foregone, since 
the exact value of T)xoii '^'^^ ensuing series development 

affords more simple results than the inaccurate one. 
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the first natural function of the "bar oscillation in pres- 
ence of a longitudinal force. 

The system f]xoj^ orthogonal and any function may 

he expanded in convergent series conformahly. (See V,l.) 

Equations (62) and (63) in conjunction with multipli- 
cation of tixqh integration of x = 0 to I, give: 



Aji sin CRa = r^- €(l-cos pr) 



dx 



(66) 



a . P / "^^on 
An cos <Pn=ira^p~ ' » " 



•n, dx 



on 



./ 'Aon 



(67) 



The addition of the squares of these terms gives the coef- 
ficients of the series expansion at 

I 



a 
1-a 



An= r-r € 



(1-cos p 



— ^ p3 ■ / ^'xon -Hx <i3C 
t) + --- sin^ pT (68) 



Pon 



2 dx 



Xon 



and the division of (66) hy (69) gives the phase displace- 
ment cpj^ at 



Pori 1 - COS p T 

tan cPn = -^^ 

" p sin p T 



(69) 



The evaluation of the integrals is hest effected with 
the exponential form of the functions Tixo^ '^^* 

Omitting the n indices, we have: 



XrtX 



/ -n^^ Ti^ dx =\'^ (l+Ho)e 



4 



Itfck- ^ Ink- e"^^ - 



1-n ^ _ l+n, 

e + — e 
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+ (l-|j,o ) e 



-(l-|j.oi)e 



1+[X 

X+\q i 



\-\q i 



"Xx. 



e 



iCX+Xo) 



e 



-(l+lJLoi)Q 



X-X^i 



Xx 



|- iXx T -i^^c 



-^0 



Taking into consideration (32), (33) expressed in ex- 
ponential form: 



XI , , -\l , 

e (1 + ij,) + e (1 - (J.) = 0 



±Xl , 
e r 1 - iJL 



and the corresponding terms for function t\j^ 

(1 - M-o) = 0 



» ° (1 + lJ.o) + e ° 



on 



iX^l 



(1 - M-o 3-)+ 



(1 + i) = 0 



(70) 
(71) 

(72) 
(73) 



The integral in the numerator of (68) reduces to 

I 

'Hxq "Hx = 2 Xq iX\j.t^XQ\x) 
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Tlxe integral in the denominator of (68) is: 

2 J 1 fl r/n, \^ a^o* / -1 \S "sXa^ 

-d-M-oi)^ i +(l+M.oO^ i e-2^^*3 + 



to : 



1 


- i 


(l+li-o) 




1 


- i 


(i~y.o) 


Cl+iXoi) 


1 


+ i 



With a view to (72) and (73) this expression reduces 
The coefficients of the series are obtained from 



■4-n = 



2 a t 



1 - a I 



€ (Xy. - Xon M-on) ( a ^^3" - 

I V '^on ~ 



vs / - cos pT)^ + sinS pT (74) 

?^on - ^ ^ V Pon 



At this time it is pointed out that all formulas de- 
veloped in this section can be forthwith generalized for the 
case of any original "bar form hy simply inserting the Z 

sign at the right-hand side of (42) to (46), (56), (59),. 
(62), (63), (66), (68), and-(74)., and giving \,A,iJ, and p 
the index n. n =1,2,3,,,. correspond to the fundamen- 
tal oscillation, first, second, ... to the upper oscilla- 
tion. However, for the reasons advanced in 71,1, the in- 
vestigation is confined to our specially selected har form. 
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Referring the shock period t to the oscillation pe*- 
riod T of the free fnndamental oscillatlont we o'btain 
with 

■ T = h T = 1) (75) 

^01 

and with due consideration to (55): 

?-2 



p T = 2 TT h ~ /l -"a (76) 
Xoi 



This enables us to compute the coefficients Aj^ for 
different a and "b from (74). (See tahle I.) The a 
values denote the magnitude of the shock load and the "b 
values the shock: period. 

The maximum h values are those at which the time 
function of (57) reaches first its maximum, iae., at which 
the maximum deflections are reached during the shock. 
This occurs according to (57) for pT = tt, or according 
to (76) when 

h = = (77) 

The shock period "b from (77) is plotted against aj; in 
figure 16, Given the coefficients for the range 0 < pT < 
TT they are known for any values of pT, inasmuch as they 
are periodic in pT. 

A 

The figures 100 cited in the last three columns 

of tahle I give the amount of the particular coefficient 
in percent of the first coefficient. The coefficients are 
seen to converge very rapidly for small a and h values 
and somewhat more slowly for high a and t values* In 
the range of 0 < a < 0,50 the first three coefficients 
suffice. The amox^nt of the third coefficient is, at the 
most, 0.3 percent of the first coefficient. In the 0,5 
< a < 0.99 range the first four coefficients are suffi- 
cient,... The fourth coefficient is, at the most, 0*2 per- 
cent, of that of the first, Por computing the deflections 
even a lower coefficient would suffice, "but for computing 
the moments the coefficients Qf. higher order grow in sig- 
nificance as their factors /^on increase with n» 
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We show lu tatle II the pha^e displacements tpji as 
computed fr^m (69). 

Taken atsolxite the maximum moment is decisive for the 
estimation of the hehavior of the tar in huckling due to 
shock load*. The 1/EJ time moments, during the actual 
shock period are, according to (56): 

& = r~-~a ^ ^ (1 - COS p t) ch X X + 

+ cos \ X + |ji Qg- Bh \ X + ^ Bin ^ a^)] ('8) 

and a fter , the actual shock period according to (59): 

--g- = 2 \on sin CPon( *-'^)+9n 3 (ch ^on 3c+°Qs \on == - 
ox n 

- sh Xon x-sin \on ^) (79) 

The maximum curvature now occurs eitheip during 

o^ after the actual shock period. It can occur during the 

actual shock only when the period lasts until the highest 
possible deflection is reached, that is, when h has at 
least reached the value given in (77)a If this happens 
the time maximum of the term denoting the period durin g 
shock is, according to (78) and (36): 



max —-^ = ----- e \^ [(1 - a)(ch \ x + u, sh \ x) + 
6x 1 - a 



+ cos K X + |JL Jl - a sin X x] (80) 
The place maximum of — :? (taken ahsolute) occurs 

when the place function 7— 8""^ that is, the hraoketed 

dx 

term in (80) reaches its maximum or minimum. The maximum 

of (1 - a) (ch \x + [J, sh A,x) is reached when x = 0, he- 

causep, < 0, Besides, since according to figure 11, 

jj, ^1 — a is a true negative fraction, positive 



cos Xx + |jL a/i - a sin A,x 
can only occur within the lengths 
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0 < X x. < 



and 



^TT<Xx<\l =4,5 approximately. 



5 
4 

Within these lengths the maximum of 



cos Xx + y, yi - a sin Xx 



likewise occurs at x = 0. The maximum of ~7~5" thus 
■ ; ' dx 

occurs at x = 0 and equals 2 - a* Figure 12 shows the 

: d Tlx 

maximum —r — = 2 - a plotted against a. The minimum of 

d? r\ — 
the place function lies in the zone of cos Xx < 0 

and sin Xx > 0, that is, ^ < "Xx < tt. It is more eas- 
ily found when assuming that the effect of the two hyper- 

d^n^ 

"bolic terms of ~ is small against the two trigonomet- 

dx^ 

rical terms when approaching the minimum. From 

• (cos X X +' \1 J\ - a sin X x) * 
dx ^ 

= "X^ (- sin X X + p, J\ - a cos X x) = 0 
the minimum is then ohtained in the neighhorhocrd of 

X X = tan"^ IJ. y 1 - a 

^^'Hx 

The computed minima of place function ~^^» taken 
absolute , are shown in figure 12. (The figures lie around 

^ = 0.65.) They are seen to "be lower than the maximum 
values for 0 < a < 0.975, so that in this range the max- 

imum occurs at x=0, For l>a> 0.975 the min- 

3x cL^^ 

imum of — r-f is decisive. For a - 0.99 its absolute 
ax 

value is only 1.8 percent lower than the maximum for the 
same a value, so the assumed absolute maximum of the mo- 
ment during , the actual shock |)eriod will be accurate enough 

->» X^ 1 1 

for all as 0.99 and b ^ =^ when made propor- 

~ X® 
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t^ohal to 
(49). to 



_3 

d T) 



= 2 



max 



a( 2__- 



a) 



Bx' 



= - 11.72' 



aA. 



max 



1 + \i \ 



or because of 



(81) 



^3 

3 ri 



The maximum of after the actual shock period 

must be defined "by trial from (79) which, filthough not un- 
dxi.ly difficult for a specific case, requires nevertheless 
a large amount of paper work for such general investiga- 
tion's as the present report, For- this reason we determine 
^bnly the upper limit value for maximum • ^ti/Bx^. 

<.:::■■: Ffj^g place-maximum of the individual Tlxon ("^9) is 

2 and occurs with x = 0. The time-maximum or mininrum of 
sin [Pon (t - t) + cpj^] is ■ ±1. As the Aj^ coefficients, 

.with .the exception of the first, are negative, the assumed 
absolute maximum of (79) cannot exceed 



3y= 



max 



3 00 2 

2 L-^i Aoi - E An A-onJ 
n=a 



(82) 



The upper maximum in (82) is almost reached, at least 
theoretically, according to thei f.ollo,wing. With due con- 
sideration to Xq-^I (see T,2,'b), (,55) gives the ratio of 
froquGncios of the (n-l)th free higher oscillation to the 
free fundamental oscillation almost exactly at 



?o.n 
Poi 



Un_+_1)_' 
" 25 



Within the period of 25 free fundamental oscillations, the 
har thus executes almost precisely: 



and 



81 1st 
169 2d 
289 3d 



free higher oscillations 
II II II 



that is, produces a beat whose period is 25 T. 

al of 12.5 T after the actual shock 

thus reaches its minimum deviation once relative 



Within the 

time interval of 12,5 T after the actual shock the value 

m 

to its upper maximum (82), In view of the great number of 
oscillation combinations during this interval and the fact 



max 
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that the frequencies increase with n, it may be assumed 



that the mazimoa of 

j 

almost or ezactly 



occurs at an instant in which 



sin [poi (t - t) - cpi] = *1 

and" sin [pon (* - t) - cpjj = + 1 are valid. 

Here it should "be pointed out that in reality certain 
discrepancies from the theoretical upper limit value given 
in (82) are nevertheless to he expected; that isi in the 
cases in which the theoretical minimum discrepancy from it 
would occur only after a greater number of fundamental os- 
cillations. For, in fact, the degree of damping is already 
appreciable and reduces the amplitudes* However, we pass 
ov.er the damping effect, as stated in the beginning. 

The cases in which the maximum ~s- and those in 

which its upper limit value a fter the actual shock are de- 
cisive, are analyzed in VI I, 2. 

b) Shock Load Greater than Eulerian Iioad 
a > 1; ^ l). 

With the initial conditions T)^ = 0 , = 0 for 

t = 0, tlie solution in this case is analogous to the sec~ 
ond Euler case (reference 1): 

T)t = r~^-- € (ch p t - 1) (83) 



where 



J OF '^^ " ^ ~ t"^ y a - 1 ' (84) 



The course of the further investigation agrees with 
that for a < 1, The equation of the elastic line for the 
actual shock period (t ^ t) is: 

= € (ch p t - 1) (85) 

a - 1 



At the end of the shock (t = t), we have: 
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- : Tl = j^-ri: ^ Cell p T - 1) Tlx (86) 

S- I?= a-^P Bhp T (87) 

whereby 'r\j^ is given in (65). 

For the condition after the shock (t = t) 
• (59) and for t = T (60) and (61) are appli,cable. 

The Ajj coefficients and the phase displacement cpj^ 
of the series (59) are obtained when equating the right- 
hand sides of (86) and (60), respectively, (87) and (61) 

with 

= —-3; e -f^ (J^M. - Kn \^on) 



V\. on 



- V^---^^) /(oil PT -if + Sh? pT (88) 

on ^ ^on ■^on 



and Pqj. ch p T - 1 

tan can = Z^—^" (89) 

p sh p T 



making 



pT=2TTh^/a-l (90) 

-^0 1 



applicable because of (87) and (84). 

The Aji coefficients have been tabulated for differ- 
ent shock loads (a) and periods (b) in table III. The 
An 

100 -r= figures -in the last three columns give the amount 

Ai 

of the respective coefficient in percent of the first coef 
f icient. 

The coefficients converge quite well, and so much the 
closer as a and b are smaller. The first four terms 
suffice throughout the range investigated. The amount of 
the fourth coefficient is, at the most, 2 percent of that 
of the first. 

Table IV gives for different a and b. 
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The time function in (85) "being aperiodic, the maxi- 
mum deflection durint^ the actual shock period is always 
reached with t = t. This instant, however, is also given 
"by (59) which decides the condition after the shock. Bear- 
ing in mind that, contrary to range a.< 1, the time func- 
tion of the speed (87) is always positive at the instant 
the shock stops and that, therefore, the deflections and 
moments continue after the end of the shock, it is readily 
seen that the maximum moment always occurs a f.te r the actu- 
al shock; consequently, the investigation may he restrict- 
ed to that of the l/EJ-fold moments of (79). 

As concerns the upper limit value of (79), the same 
arguments and the same formulas are applicahle as for 
a < 1, that is, the upper limit value of (79) is given "by 
(82). It occurs at' x = 0 in spite of the fact that - 
as seen in f igiire 12 - for a > 1 at the instant of in- 
cipient free oscillations the maximum value of the moment 
for X = 0 is smaller than its assumed ahsolute minimum 
value occurring at around x = 0.65 1. 



The ohvious criterion of the mas'^uitude of moments in 
huckling due to shock loads is the moment of the same har 
in static "buckling stress. For this reason we shall at 
first, attempt to resolve the differential equation for 
static "buckling of a "bar of proscrihod form. Thus we of- 
foct the calculation for tho "bar conformahly to (46), for 
-example, although it could bo given direct for this par- 
ticular case, as shown elsewhere. 

The differential equation of static "buckling is af- 
forded from (8) when omitting the term with a^y/Bt^, al- 
though figure 13 reveals it even more elementary. It is 



VII. RESULTS AlID INTEEPEETAT lOH 



1* Comparison with Static Case 



M = P (ti + y - y^ax) 



A (I 



x) 



and consequently, 



a^ri _^ P 1 



EJ ' EJ 



[A (I - x) 



P (y 



ymax) ] 



(91) 
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Ihe homogeneous equation, 

resolyes to 

1^ = 01 cos CO X + C3 sin co a: (93) 
Putting (93) in (92) gives: 

and consequently, 

0)2 =1? - \* (94) 

"because of (16) and P = - X, Equation (91) is resolved 
"by variation of the constants in (93), that is, from 



= Ci(x) cos w X + Cs(x) sin CO x (95) 



where 



, ^ f(x) Ug dx 
Ci(3c) = / —7 —7 + ^1 

O 111' - Tig' 111 

X f(x) u, dx 
c(x) = / — -l-i—i— + kg 

O Us' Ui - Ui« Us 

_ f (^) = - - ^ ^7 - y^ax)] 

U^ = COS (0 X 

Us = sin <o X 

and kj^ and kg are constant values, The terms for y 
and ymax are taken from (46) and (47), After appropriate 

reductions and calculation of ' and Cs , it gives 
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4 .- 1 



^ = -.(1-3: - Icoscox. -(jj sia, Oi)-x) - 



- y + yjnax^-"- " cos CO x) - c rcli A x+-^cos \ x + 

— 3 

+IJ, (^sh X X + ^ sin X. - oo x + cos o) x^ 

+ kx cos CO X + kg sin (o x (96) 

Tlie ■boundary condition r\ = 0 for x = 0 gives 

ki = 0; those of ^ =: x = O and ri = 0 for x = 0 and 

dx 

A 

X = I give two equations for - and kg, thus yielding 
with allowance for (32) and (33): 

2 A' ' . .. 
ka = - = 

and 



A 

The introduction of .the values of kj_ , kg, and — in 
(96) with consideration to (36) and (94) gives: 

= e — [ch X x-cos ?v x+ijl (sh X x f ^ sin A, x)] (97) 
1 - a • ^ ■ 

This term approaches oo as a approaches the value 1. 
Consequently, the Eulerian "buckling load occurs at a = 

-gr =1, as already explained in VI, 3. Since then, hi = 

EJA. : -.• 

4.5 (ahout) according to figure 4, we ohtain for the Eu- 
lorian load the known valuo 

Ptj = 4.5^ li '5 2 TT^ M (98) 

independent from the original "bar form. 

The l/EJ-fold moments aro chtair-ed in first approx- 
imation from (97) at 
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4r ch X X + cos X X + 



d?" 1 - a 

+ p, sh \. X + sin \ x)j (99) 

Comparison with (80) reveals them to he exactly half 
as high as the time maximum of the l/EJ~fold moments due 
to shock load "by a time-constant force, for an effective 
period equal to half the period of the oscillations accom- 
panying this load.* 

This result could equally have heen obtained direct: 
The solution (43) of (8) represents at the same time the 

resolution of (8) without the p P term, when equat- 

at* 

ing Tj-t = k = constant* Then (44) "becomes: 

C (E J + X) "X® k + X G \® ] Tlx = 0 . 



hence , 



€^_X 

B + X 



1 - a 



Or, hy the same argument, so long as the shock force 
is constant in time and lower than the Eulerian load, the 
oscillations are harmonic* If the har form is such as to 
insure only one definite oscillation - in the present case, 
the fundamental - the mean position of the oscillation cor- 
responds to the equilihrium position assumed hy the bar 
under static load. The mean position is reached after a 
one fourth period} the position of potential total energy 
and consequently of maximum deflections and moments is 
reached after one half period. The oscillation helng har- 
monic, tho maximum deflections and moments are exactly 
twice as groat as the corresponding magnitudes in the moan 
position. 

In tho general case, however, with the har form such 
as to produce a motion consisting of tho superposition of 
sevoral normal oscillations, the motion of tho individual 



♦That stresses due to "sudden" loads can "become twice as • 
high as in static loading, is a well-known fact, (Ponce- 
let, Introduction a la Mecanique Industrielle, physique et 
experimental e, 1830.) 
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"bar particles is not synchronous, nor is there any condi- 
tion in whici. the total energy has potential form. In this 
case then there is no siraple relationship "between the de- 
flections and; moments prodticed during the period of load- 
iing due to static and to "sudden" loading.* The deflec- 
tions and moments for the static case are then ohtainod 
from the corresponding magnitudes of the static case hy 

s.v.'bstltuting . --- — £n for the time function t\^j^* or; hy 
direct treatment of the static case. 

These considerations referred to cases of tension 
lower than the Eulorian load (a < l)# For tensions great- 
er than tho Eulorian load (a < l) the approximate equa- 
tion of tho clastic lino fails, as known (curvature = 
d y/dx^). For duch cases the exact equation should he 
used, i.e., write: , • . . 

d ri 

---T *k-^ = h ^^'^ -i"<y-y^:c)3 . (100) 

instead of (91). 

.The resolution. of (lOO) is very difficult. Even a 
second, approximation as attempted hy Trefftz ( reference:- 5) 
seems intractahl©. in view of the complicated interference 
. function of ( 100) . Added to this is the fact that when ef« 
footing the dynamic analysisj we also substituted its 
first approximative value d ri/ dx^ for the curVature, A 
comparison of the results of the dynamic study obtained in 
first approximation with those of the static study in sec- 
ond approximation would, however, not he logical* 



•These ■ explanations aro, according to Section 128 of lovo- 
Timpe's volume on ela.sticity. It states, among other 
things that, in the general case "The distortion remains 
less than twice the distortion existing in oquilihrium po- 
sition." This docs r.ot seem strictly valid. It is well 
posisihio- that at isolated points the distdrtion may he 
greater than tho twofold distortion existing in static equi- 
lihriua position,. cve:i if tho energy of tho whole har, 
whi-Qh is governed hy an integral term along the bar, is 
not pi3.rely potential at any instant.. . (Compare VII ,2 of 
this report, which contains similar deductions for tho 
free oscillations following the actual shock.) 
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. Under these circumstances we content ourselves with 
iCQmpariBig the static to the dynamic case for a < 1, while 
for evaluating the results of the dynamic, case for a> 1, 
we only adduce that of tho dynamic case for a< 1 in 
comparison* 

2, Numorical Costparison, Doduct.ionst and Sxamples 

.On tho hasis of the foregoing, the function for the 
'moment in tho static case for a < 1 is, according to (81) 



dx 



a \ 

max l + jx^-l 



- 5.86 —7-— —r ^o (101) 



The ratio of the maximum moments (taken absolute) in 
the static and dynamic case is theni according to (Sl) and 
(101) for a < 1:. 

c = 2 (103) 

a) Provided the maximum moment occurs during the. ac- 
tual shook, and 

h) Provided the upper limit of c is according to 
(82) and (101): 

. = _ 0^*|a_t^^_kl [i, - X ^ ^"onl (103) 
Eq a A n=s 

?rom (102) and (103) it is seen that in our case the 
assumed absolute ratio of the decisive moments in the dy- 
namic and static case is independent of the amount of orig- 
inal ecGLentricity Cq (A^/Co does not contain €^)« 

The evaluation of (103) affords the figures 14 to 20. 
The particular talhlSB have been omitted. Suffice it to say 
that taking' into account the terms of the series expansion 
applied in table I, the amount of tho term of maximum or- 
der is, at the most, 2 percent of that of tho first term, 

Figure 17 gives tho c values as ordinates, the b 
Values as abscissas in logarithmic scale. The a and a.-^ 
values (in parenthesis) are the parameters. The range of 
each carve was chosen so that the maximum shock period (b) 
corresponds to the maximam deflection during the actual 
shock period. (bjnax tlien corresponds to (77).) 
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• '■ ■ ' • , ■ ; : 

According to fignr© 14 the stresses in short shock pe- 
ri-6d^" are lower than for an identical •;Btat?,_o.. stress. But 
^'i^fc-' longer shock periods the static stresses can "be exc^ed-r 
elL 'considerably. 

The maximum c values occur at shock periods corre** 
spending to (77). They are, on the whole, greater than the 
maximum c value of (102) for the same shock period. 
Therefore, the maximum moments are even then only to he ex- 
pected after the actual shock period, when the greatest de- 
' -flection possible in this period is reached during the actu- 
al shock period. 

This somewhat confusing statement needs further expla- 
nation. If the actual shock lasts just long enough to en- 
sure the greatest deflection so that, as in our case only 
one definite normal oscillation occurs, then the total en- 
ergy at the end of the shock is existent in the form of 
deformation energy. During the ffee oscillations* follow- 
ing the shock there is, therefore, no instant in which the 
energy of deformation is greater than at the instant of 
shock termination. Generally the first is even consistent- 
ly lower than the latter for, aside from the type of sup- 
port according to the second Eulor case, there is no inte- 
gral ratio between fundamental and upper oscillation fre- 
quency, hence no instant during the free oscillations in 
which the total energy is potential. 

The result, ©.ccording to which the upper boundary 
value of the maxim am moments after the shock is greater 
than during the shook at the moment of maximum deflection, 
is., however, readily understood when bearing in mind that 
the highest moment in a certain period does not necessari- 
ly have to occur at the instant when the total energy is 
energy of deformation. Hather, there may be instants in 
which the superposition of oscillations may even induce 
still higher moments, even though the concurrently exist- 
ing deformation enyrgy of the whole bar - which, as is 
known, is proporb ional to the integral of the square of 
the moment - is lower than the doforraation energy at the 
end of the shockr 

Figure 15, which is simply the modified figure 14, 
reveals that in the range of higher shock loads the dynam- 
ic stresses pertaining to the same shock periods increase 
more slowly with increased shock load than the correspond- 
ing static stresses (and vice versa). 

♦Excepting the second Euier case, several free oscillations 
occur as a rule after the shock even when only one defi- 
nite oscillation existed during the shock. 
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Contrariwise, figure 16 shows that the greatest pos*" 
sihle c values increase with the shock loads wherehy, 
however, the shock periods (h) for these values increase 
equally. In our case, when the shock period is long 
enough, the dynamic load may reach around 3.4 times that 
of the static stress, whereas with vanishing shock load 
.the stress ia at the most, twice as high as the static 
stress's 

Figure 17 shows the shock periods (h) against shock 
load (ag) for shock stress equal to static stress. 
These shock period's are seen to increase slowly at first, 
then faster, and to approach infinity in the neighhorhood 
of the Eulerian load. 

Figure 18 shows c against shock period for as = 
0*770, together with periods exceeding those of (777. 
After "becoming maximum, c at first drops with increasing 
shock period hut can, of course, not drop below 2, since 
this figure had already heen reached once at the instant 
of maximum deflection during, the shock. Thus, starting 
from a stated shock period c = 2 is valid. At still 
longer shock periods, c then increases again, etc. At 
other a£ value.s, the trend of the curve is quite similar. 

To render the interparetation of the dynamic test data 
possible in the a > 1 range, in view of the absence of 
comparative base of static case, we plotted the dimension- 
less qtiantities: 



2 



l_ ojn 



[Ai \oi - 2 An Kojil 



max 



(soG (82)), which are proportional to the upper boundary 
values of the moments for the wholo investigated range of 
a and against b in figure 19, with a and as as 

parameters, abscissas and ordinates being given in loga- 
rithmic scale to ensure better siirvey. Figure 20 gives 
the same results, but b serves as parameter and the a-^ 
values form the abscissas. 

According to figure 22, there is no maximum moment 
when a ^ 1; the moments continue to increase with the 
shock period, while for a < 1 the moments reach a maxi- 
mum which may not be exceeded during any shock period. 



The Eulorian buckling load (as = l) thus represents 
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a remarka'ble iDOtiiidary for the case of shock load. But, 
while denoting a stahility limit for the static case, "be- 
yond which the har fails,* it may he safely exceeded in 
dynamic load provided the period of loading is short enough 
and the crushing limit of the material is not exceeded 
therehy. To illustrate: If in our particular case the 
shock period is 0*6 times the period of the free fundamen- 
tal oscillation of the har (h = 0*6), the maximum stress 
\xpon exceeding the Eulor load hy ahout 17 percent (ajj = 
1,167) is only ahout as high as the highest possihlo 
stress under a shock load, which is ahout 14 percent (au 
= 0.864) below the Suler load whereby, however, the lat- 
ter is reached only with a longer shock period (h = 1,2), 

To gain an insight into the order of magnitudes of 
the shock periods T (equivalent to the range of h val- 
ues in this report) and the allowahle maximum "relative 

original eccentricitios" (5,8iS -—)** several examples 

have "boon computed in table T. 

Prom (75) and (54) follows: 



Besides, since our investigation applies only to con- 
ditions in which t?ae stresses must remain below the pro- 
portionality limit, it follows that 



which may also be expressed as 



and the highest permissible "relative original eccentric- 
ity" is: 



*The approximate differential equation of the problem is 
based upon these arguments. (See VII, 1.) 

**Soo equation (48), 




(104) 




P 
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5*86 "5" 5*86 —p 

I - EJ 1 



The -r- 



8 Tl 
8?" 



(105) 



max 



max 



values are given in figare 19, 



TaTjle V gives (104) and (105) for a normal steel tuTae 
of 48/l mm diameter with a proportionality limit of Op = 
2,250 kg/cm^, for a high-tensile 24/l mm diameter steel 
tube with <7p = 4,200 kg/ cm^ proportionality limit and a 
48/l mm diameter duralumin tuhe with CJp = 1,400 kg/ cm^ , 
each for a load stage 23 percent helow and 16*7 percent 
ahove the Eulerian load.* Two har lengths were considered 
for the lowor stage; the shortest corresponds to the mini- 
mum slendernoss 



I = 4,5 




P 



to which the Eulerian load is yet applicable (compare (98)), 
the longest one to |- = 200. At the upper stage only 

~ = 200 was considered, since the concept: boundary slen- 
dernoss conformal to Euler, is without meaning. 



It is seen that the permissible eccentricities de- 
crease as the shock periods and loads increase. A crite- 
rion for the evaluation of the computed eccentricities is 
Muller-Breslau' s recommended = 0.005. 

The order of magnitude of the shock periods in the 
analyzed examples is 0.001 to 0.07 second. 

In reality it should be possible to exceed the eccen- 
tricities of (105) without danger, for tho elongation and 
elasticity limit decisive for the dynamic loads is proba- 
bly higher than the corresponding limits for static loads 
(reference 8), The execution of systematic studies for 
the determination of the limit of dynamic elongation and 



*The Op and B values are taken from A. Rechtlich's re- 
port^, published in the 1931 DVL Yearbook, p. 379; and 
Schroder* s report as published in the 1928 DVL Yearbook, 
p. 216. 
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elasticity of materials used in bridge and airplane design 
should "be vrell worth while. 

Likewise it is very desira'ole- to checlc the results of 
this study hy actual test. The writer hopes to te ahle 
to do this in the near future. 

5. Range of Validity of the Approximate 
..■ . - , Differential Equation 

A check on whether and to what extent the use of the 
equation in which f ) - was disregarded, is admissible. 
To this end we attempt to determine whether the (§"j^) 
Values are actually small against 1, (See IV, 5.) 

The check may he confined to the maximum deflections 
after the actual shock, since, as shown in VII, 2, they are 
nearly always greater than those during the shock. 

Equation (59) gives: 
|5 = g ^n \)n sin (po^ t + cp) ( sh X^^ x 

+ sin' X - ch \on x + cos x) (106) 

The individual place functions in this sum are maximum or 
minimum for 



d 3 

--p-- = (ch Xon + cos Xon x 

- sh X - sin Xon x) = 0 



dx^ 



e" ° = sin Xq X ~ cos Xq x (107) 

when omitting the indices n and considering 

, ^ ~ Xo X 

ch Aq X - sh A.0 s = o 

Equation (107) is resolved with 

Xq- X: = 1.036,' 3,91, ~ S TT, ~ i~ TT. ^ 5:^^ TT 

4 4 4 

c 

For the fundamental oscillation Xq I - -g n is appli- 
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ca"ble. Consequently, ^ Tjecoiaes marinnim either fop 

ac ^ 4Ji_l,.g36 ^ 0,2 63 
I 5 TT 

or 

X = 4_X_3a91 ^ 0,995 
I 5 TT 

The evaluation gives: 

?!^= 1.019 for S = 0.263 

■ . d X t- 

^ = - 1.434 for S = 0.995 
d X I* 

Hence ^ = 0.995 is decisive for the fundamental oscilla- 
tion. We substitute here ^ = 1.00 for which, likewiso, 

— "Q-S = approximately - 1.434 is valid, 
d X 

But for 5 = 1,00, the place functions of all higher 
X 

oscillations also reach a maximum or minimum, since the 
higher solutions of (107) agree with t3ie higher proper 
values of our differential equation. Specifically, It is 



max 



^Tlxn 
dx 



~ sin ^^tl ^ + cos TT=8 sin J = 1.414 



— X X 

since e ° = ch 3C - sh >o 3C 

hocomes nogligihly small for higher values of Xq x. 

By the same arguments as used in VI, 3, a, we ohtain as 
upper limit value for max i^ftn/^xj : 

list " ^•434 Ai \oi - 1.414 Aa \on (108) 

IS'O i X 
T~ can he computed from (108) by means of tables 
o A Co 

I and III, which contain the £~ values. For -r~ the 

permissible eccentricities (see VII, 2) must be included. 
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This calculation was carried th.rougli for tho exanrplos 
of tablo T and compiled -in ta"blo VIn . 

Co 

Limiting the eccentricity so as to keep the stress- 

es during the shock Td el ow the proportionality limit, the 
i~\ values for a stated "bar and stated shock load scarce- 



X ! 

ly depend on tho shock period; furthermore*, tho 



val~ 



ues alsove the Eulerian limit. are even-smaller than those 
helow it, 

a 

In the examples the maximum, is (^^^ "~ 0«082^ = 0.007 

for a high-tensile steel tube (24 X 1 mm» I = 163 cm, 
p-g = 750 kg) at a = 0.77 and^ "b = 0.10, that is, for a 
comparatively small shock load and ^hock period. Even 

this max {^^-^J "be summarily ' disregarded against 1, 

The use of the approximate differential equation is 
herewith proved as justified, 

APPElIDiX 



Derivation of the interference functions for the sec- 
ond and third Euler'.casb. 

For tho pin~ondod support as shown in figure 21, we 
have; for x ~ I the differential equation of the elas- 
tic line, hecauso of 

, . M = A X + P n = -'^-v--^ X + P -ri 

at ^2 

E J + P n 



+ p ^ = . 
which resolved, Ijecomes: 

I - g 

't) = A sin w X + B cos w X - —-^ — x 

where 

The limit condition ri = 0 for x = 0 gives B = 0. Con- 
sequently, 
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: r\ ^ A aia o) x - x- '• • (109) 

Pi. 



Tor .3? =.t .tlia differential equation ia.: 

^ ^ + ^- = - ~f ~^ 3= + (-3= - 

and. its resolution Is: . 

Ti = A sin CO X + S COB 0) X - ^u-i^-S | 

The limit condition = 0 for x = I gives: 
B .= - A tan CO I, 

consequently, 

T) = A (sin (0 X - tan (o I,, cos oo x) - — i (llO) 

The determination of constants A and A follows from 
the continuity conditions at x = | .. Per this pcint the 
n values of (109) and.(llO)', as well as their first de- 
rivatives, must "be mutually eqUal. Equating (109) and. 
(110) for xw ^ gives: 

A = A (1 - tan CO I cot CO |) (ill) 

Equating the first derivative gives: 

A cos CO I - A (cos CO I + tan co I sin <^ i) = ^ (112) 

Then (111) and (112) yield: 



A = - 



si n CO i 

F CO tan (O I 



A = sin CO il - i) 
F CO sin CO I 

and (109) and (llO) hecome: 

Kl (x.O = |i^JP4L_-_il Bin CO X - 1-t-^x for x< i 

Fcoamcot PI " 

Kj.(x.|) = |i5^^_LL_r_£l sin CO I - ^^i^i " t 

Pcosincol Pi 
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where ve use r| in place of Z{xti)* As seen* the core 
ftinctions in z and i .actually are symmetrical. 

At the expense of much more paper •vrork, although fun- 
damentally with the same degree of accuracy as in the ca;se 
of figure 21, tho core functions may "b© obtained for other 
support conditions* Tor the tar clamped at one end and 
hinged at the other, for example (see fig, 22), the un- 
known support pressure 4. occurs in addition to the four 
constants A, S, A, and S« The conditions for resolving 
these five quantities are: 

-n^ = 0 for X =; 0 

= 0 " X = I 

'r 

-n^ = " X = I 

Tly' = 0 " X = I 

■■ 'Hr' = 'Hi " ^ =■- ^ 

wherein suhscripts - I and r denote the parts of the 
elastic line on this loft- and right-hand side of x = |. 
The calculation gives the following symmetrical core func- 
tions: 



For X = I 

1 I sin uj x+x sin uj|- ujx | cos lu I - |^sin(ju (|- Z) +U)Z cos tu^;^ -!)]] sin cu x x 

= - sin lar- 0,1 cos 0)1 '■ P 

for X ^ I 

1 X sin u)|+| sin ujx-oux | cos ujZ - - [sinu;(x-l) +u;Z cos tu(*- Z) Jsin cu| | 

Kr(x,|) = ; — — = - 

P sin cuZ - cuZ cos tJuZ . P 
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TABLE I* Coefficients An of Series (59) for a < 1 





a ' 


t 










100 ^ 










i 










Ax 








n=l 


n=2 


n=3 


n=4 


n=2 


n=3 
















p e r c e 


XL t 


0. 


25' 


0.10 


0.159 


—0 . 00052 


—0 .OOOOR 




0,33 


0,05 








0.30 


0.453 


-0.00282 


-0 ,00054 




0.62 


0.12 








0.556 


0.632 


-0.00488 


-0,00095 




0,77 


0.15 




0. 


50 


0.10 






— 0 OOfT^fi 

— w aWwOO 




0.75 


0.10 








0 .25 


V # O 9 u 


• U JL J> W 


«o om RR 




1.26 


0.21 








0 . 50 










1.72 


0.30 








0.651 


1.863 


-0.0332 


-0,00572 




1.78 


0.31 




0.75 


0 . 10 


0, 571 


-0 .00746 


-0 .00096i 


-0 .00025 


1.31 


0.17 


0 .04, 






0 . 30 


1 .962 


-0 .04595 


-0-00724 


-0 .00199 


2.34 


0,37 


0.10 






0 .50 


3. 59 5 


-0 . 1028 


-0 .0165S 


— 0 .00460 


2.86 


0.4G 








0 .859 


5.232 


-0 .1622 


—0 .0263fi 


—0 .0073'', 


3.11 


_0^50 


O 14 

\J . X'X 


0.84 


0 .10 


0.680 


-0.0103 


-0.00132 


-0.00028 


X . ox 


U • 


0.04 










— A A A »5 A 
—0 .04^0 


— A AA*90K 


.-A AAl Q '7 






A T A 






0 . 50 


4.755 


-0 .1533 


-0.0246 


—0 .00 690 


3.2? 




0.15 






0.75 


7. 565 


—Qt .2 66 


-0 .0431 


-0 .0120 


3.48 


0.57 


0.16 






1.040 


8 .972 


— 0 .32? 


-0 .0520 


—0 .0146 


3. 58 


0 . 58 


0.16 


0.91 


0.10 


0.771 


-0,0137 


-0.00167 


-0,00042 


1.78 


0.22 


0.05 






0.25 


2,245 


-0.064 


-0.00939 


-0.00247 


2. 84 


0,42 


0,11 






0.50 


5.94 


-0.222 


-0.0340 


-0.00902 


3,74 


0,57 


0,15 






0.75 


10 . 50 


-0 . 426 


-0.0657 


-0,01745 


4,06 


0.63 


0.17 






Jl aUvJ 


X It* Do 


— A A 


A AO AA 


..A A OKA A 


4.18 


0. 64 








J. . o4 




_A C 
— U , f J. D 


—A 1 T A C 


_A AOQAA 
— U .Uos^bU 


4.24 


0 .65 






95 


0.10 


0.819 


-0.0159 


-0,00189 


-0.00047 


1.94 


0.23 


0.06 






0.25 


2.425 


-0,0748 


-0.010 62 


-0,00280 


3,08 


0.44 


0.12 






0.50 


6.695 


-0.2695 


-0.0399 


-0.01060 


4.02 


0.60 


0.15 






0.75 


12.64 


-0.5515 


-0.0823 


-0.02195 


4.36 


0.65 


0.17 






1.00 


19,22 


-0,867 


-0,129 6 


-0,0345 


4.51 


0.67 


0.18 






1.50 


29.43 


-1.355 


-0.2027 


-0.0541 


4,60 


0.69 


0.18 






1.764 


31.03 


-1,432 


-0.2144 


-0.0 572 


4,61 


0.69 


0.19 


o7 


99' 


' 0.10 


0.877" 


-0.0185 


-0,0022 


-0 ,00043 


2.11 


0.25' 


"0.05 






0.50 


7.67 


-0.3333 


-0.0485 


-0.0105 


4.34 


0.63 


0.14 






1.00 


26.00 


-1.257 


-0,185 


-0.0403 


4.84 


0.71 


0.16 






1.50 


53.05 


-2,622 


-0.388 


-0,0845 


4,95 


0.73 


0.16 






2.00 


84.50 


-4.21 


-0.623 


-0.1358 


4.98 


0.74 


0.16 






2.50 


115.20 


-5.76 


-0.852 


-0.1858 


5.00 


0,74 


0.16 






3.00 


139.9 


-7,01 


-1.038 


-0.2260 


5.02 


0.74 


0.16 






3.50 


154.9 


-7.76 


-1.149 


-0.2510 


5.02 


0.74 


0.16 






3.845 


158.0 


-7.93 


-1.173 

i 

L . 


-0.2560 


5.02 


0.74 


0«16 
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TABLE II. Phase SiBplacements cpxx of Series (59) 

for a < 1 



■.. . .a 


t 




n=l 


11=2 


n=3 


11=4 


0.25 


0.10 
0.30 
0.556 


0 , 3131 
0.8988 

Tt/2 


0 , olO 
1.330 

tt/2 


1.454 

tt/2 




0*50 


0.10 
0.25 
0.50 
0.651 


0.3097 
0.7242 
J. . c uo o 

Tr/2 


0.783 
1.224 

1 AAA 

tt/2 


1,122 
1.398 

tt/2 




0.75 


0.10 
0.30 
0.50 
0.859 


0,3074 
9,8969 
1.1489 

Tr/2 


0,799 
1,176 
1.433 

tt/2 ' 


1.134 
1.430 
1. 504 

7t/2 


1,304 
1.489 
l;532 

tt/2 


0.84 


0.10 
0.25 
0.50 
0.75 
1.038 


0.3127 
0,7004 • 

1.1111 
1.3538 
tt/2 


0.809 
1.220 
1.419 
1.505 

7t/2 


1.410 
1.414 
1.50 5 
1.543 

Tr/2 


1.309 
1.469 
1,528 
1.553 

Tr/2 


0.91 


0.1 

0.25 

0.50 

0,75 

1,00 

1.34 


0,3589 
0.6820 
1.0597 
1.2716 
1.4153 
tt/2 


0.796 
1.208 
1,400 
1.47 6 
1.523 
Tr/2 


1,131 
1.391 
1,488 
1,525 
1,598 
Tr/2 


1,303 
1,465 
1.522 
1,544 
1,557 

1,302 
1,463 
1,519 
1,540 
1.551 
1.565 
tt/2 


0.95 


0.10 
0.25 
0.50 
0.75 
1.00 
1.50 
1,764 


0. 3043 
0.6729 
1,0339 
1,2255 
1.3448 
1.5'b30 
tt/2 


0,794 
1,202 
1.389 
1.447 
1,500 
1,550 
7t/2 


1,131 
1,388 
1,483 
1,518 
1.536 
1.561 
tt/2 


0.99 


0.10 
0.50 
1.00 
1.50 
2.00 
2.50 
3.00 
3.50 
3.845 


0.3020 
1,0100 
1,2787 
1.3878 
1.4493 
1.4910 
1.5290 
1,552 
Tr/2 


0,790 
1,379 
1,478 
1.514 
1,533 
1,546 
1.554 
1.565 
tt/2 


1.112 
1,474 
1,525 
1.542 
1.552 
1.558 
1.563 
1.568 
Tr/2 


1,300 
1,516 
1.595 
1.555 
1.560 
1.564 
1.566 
1.569 
Tr/2 
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a?ABLE III. Coefficients An of Series (59) for a > 1 







■■ ■- .-. ■ 






















100 - — 




a 












Ai 








n=l 


n=2 


n=3 


n=4 


n=2 


n=3 


n=4 















P e 


r c i 

: 


9 n t 




1.01 


0.10 


0.911 


-0.0201 


-0,00230 


-0.00028 


2.21 


0.25 


0.03 




0.25 


2;768 


-0.0969 


-0.01338 


-0.00173 


0 . OU 


f\ Ad 


0,06 




0.50 


8.20 


-0.371 


-0.0532 


-0 .0070 


A f^O 


U , DO 


0.09 




0.75 




-0.843 


-0,1219 


— A Ate 


4. ( 0 


U . r X 


0 .0 9 




1«00 


OO • 4; 


— X . OoU 


—0 , eieia\J 


A A 0 Q 


n't! 




A T A 




X • DO 


r JL • D 


-3.688 


-0.536 


A A r> A 


5.16 


0,75 


A T A 




<o •UU 


Loo • f 


-7.21 


-1,048 




5.20 


0,76 


A . T A 


1*05 


0. 10 


0.974 


-0.0238 


-0,0026 


-0 .00064 


2.44 


0,27 


A A r> 

0.07 




0,25 


3.025 


-0.1163 


-0.0155 


-0.00393 


3.85 


0,51 


0.13 




0.50 


9,44 


-0.466 


-0.0645 


-0.01652 


4.94 


0,58 


0.18 




0.75 


21.47 


-1.140 


-0,1593 


-0.0408 


5.31 


0.74 


0.19 




1.00 


41.9 


-2.285 


-0,3205 


-0.0823 


5.46 


0,77 


0.20 


1.09 


0,10 


1.072 


-0.0293 


-0.00314' 


-0.00075 


2.73 


0,29 


0.07 




0.25 


3.414 


-0.1448 


-0,01858 


-0.0047 


4.24 


0,54 


0.14 




0.50 


11.36 


-0.617 


-0.0818 


-0.0209 


5.44 


0.72 


0.18 




0.75 


28.20 


-1. 629 


-0,2195 


-0.0561 


5.78 


0.78 


0.20 




1.00 


61.70 


-3. 660 


-0.4930 


-0,1262 


5.93 


0.80 


0.20 



TAB L E IV . Phase Displacemen ts cp^ of Spries (59 ) f or a > 1 



a 


t 








n 








n=l 


n=2 


n-3 


11=4 


1.01 


0.10 


0,3010 




0,788 


1,126 


1.307 




0.25 


0.5630 




1.349 


1.427 


1.464 




0.50 


0.9949 




1.452 


1.476 


1.516 




0.75 


1,1573 




1.490 


1.506 


1.534 




1.00 


1.2457 




1.609 


1.521 


1.542 




1.50 


1.3360 




1.527 


1.536 


1.551 




2,00 


1.3788 




1.535 


1,542 


1.554 


1.05 


0.10 


0.3024 




0.802 


1.127 


1.300 




0.25 


0.6571 




1,198 


1,381 


1.459 




0.50 


0.9703 




1.367 


1.470 


1.511 




0.75 


,1.1117 




1.423 


1.498 


1.538 




1.00 


1.1823 




1,448 


1,510 


1.535 


1.09 


0.10 


0.3024 




0.786 


1.130 


1.300 




0,25 


0.6473 




1.184 


1.378 


1.457 




0.50 


0 .9413 




1,350 


li464 


1.508 




0.75 


1.0630 




1,401 


1.489 


1.523 




1.00 


1.1164 




1.422 


1.499 


1.528 
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TABLE v.. Illustrative Examples of Highest Permissiljle Rela- 
tive Original Eccentricities 7jaa.x/^ ^'^^ Different 
Tutes Under Different Shock Loads P and Different 
Shock Periods t . 



Steel tulje 48/1 mra-l^ = kg/cn? 

' ia-o = 2250 II 





h 


f ^ = 134 
P]i}=3330 kg -j 1 

= 223 cm 


J i- = 200 
P]jl=1500 kg-(^ i 

I = 332 cm 


P 


T 


I 


P 


T 


I 








S 




kg 


S 




0.770 


0.10 
0.25 
0,50 

0,75 
1,04 


2560 


0,0024 
0,0061 
0,0122 

0,0183 
0,0254 


0,0126 
0.0042 
0,0016 

0.0010 
0.0009 


1160 


0.0054 
0.0136 
0.0271 

0.040 6 
0.0564 


0,0533 
0,0179 
0.0070 

0.0043, 
0.0037 


1.167 


0,10 

0.25 
0.50 

0.75 
1,00 








1755 


0.0054 
0.0136 
0.0271 

0.0406 
0.0542 


0.0234 
0.0069 
0.0020 

0.0008 
0.0003 



(mm X .03937 = in.) (cm X ,3937 = in.) 



(kg X 2.20462 = Ih.) (kg/cm^ X 14.2235 = Ih./sq.in. 
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TABLE T (continued) 



Steel tnte 24/l mm = ,2X10« kg/ cm^ 



Up 



as 


t 


Pjj=3050 kg 1 i " 

= 80 cm 


, i. = 200 
Pe=750 kg il 

= 163 cm 


P 


T 


7max 

■ ■ I 


P 


T 


ymax 
I 






kg 


8 




kg 


8 




0,770 


0.10 
0.25 
0,50 

0,75 
1.04 


2345 


0.0006 
0,0016 
0.0032 

0.0048 
0.0067 


0.0167 
0.0056 
0.0022 

0.0014 
0.0011 


565 


0.0027 
0.0067 
0.0133 

0.0200 
0.0277 


0.1214 
0.0409 
0.0160 

C.OlOO 
0.0084 


1.167 


0.10 
0.25 
0.50 

0.75 
1.00 








855 


0.0027 
0.0067 
0.0133 

0.0200 
0.0266 


0.0649 
0.0191 
0.0055 

0,0022 
0.0010 



l»tiralamin tute 48/l mm -1^^ _ 



0.75X10^ kg/cn? 
1400 " 



as 


■b 


. i - 

Pe=2060 kg \ i ' 

■•If r 


= 104 
: 173 cm 


Pe=560 kg 1 i 


200 

332 cm 






P 


T 


ymax 

~T" 


P 


T 


ymax 
I 






kg 


S 






S 






0.10 




0.0014 


0.0165 




0.0053 


0.10 68 


0.770 


0.25 
0.50 


1585 


0.0036 
0.0072 


0.0056 
0.0022 


435 


0.0132 
0.0264 


0.0360 
0.0141 




0.75 




0.0108 


0.0014 




0.0396 


0.0087 




1.04 




0.0149 


0.0011 




0.0549 


0.0074 




0.10 










0,0053 


0.0561 




0.25 










0.0132 


0.0165 


1.167 


0.50 

0.75 
1,00 








655 


0.0264 

0.0396 
0.0530 


0.0047 

0.0019 
0.0009 
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TABLE VI. Upper Limit Values of 
the Examples of laTjle V 



a 



as 



OH 



for 



Upper limit of 1 c>n/5xl for 



Steel tulje 
diameter 
48X1 mm 
Pj,=1500 kg 

I ~ 232 cm 

^ = 200 



Steel txCbe 
diameter 
24X1 mm 
Pj.=750 kg 

I =163 cm 

\ =200 



Duralumin tuTae 
diameter 
48X1 mm 
Pjj = 560 kg 

I = 332 cm 

^ = 200 



Eccentricity 

to talile V 



max 

"V 



- according 



0.84 



0.770 



0.10 
0.25 
0,60 
0.75 
1,04 



0,036 
0,034 
0.034 
0,034 
0.034 



0.082 
0*080 
0.079 
0.079 
0.075 



0.072 
0.070 
0.069 
0.068 
0.068 



1.09 



1.167 



0.10 
0.25 
0.50 
0,75 
1,00 



0.026 
0.025 
0.024 
0.024 
0.020 



0.070 
0.069 
0.067 
0,0G7 
0.067 



O.X361 
0.059 
0.058 
0.058 
0.058 



Translation "by J. Vanier, 
national Advisory Committee 
for Aeronautics. 
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Figs, la, lb, 2, 3 




Figure la.- 
The "bar 
hinged 
I ] at "both 
ends 
under 
tension P.. 



Ax 



a 



As 
A(Ti+y) 



Figure Ih.- 
Element of bar 
of figure la. 



40r 



30 



10 



20 



10- 



100 Z 200 



Figure 2. -p (boo fig. la) versus 
e 

length of bar. 




f(a,^) 



f(X) =tanM and ^(a-'^)=.,^2^ th(X ZVi-a ). (parameter a), 



f for a =1,05 
f for a =1.00 
(Siilerian load) 
f for a =0.95 
f for a =0.75 
f for a =0.5 
f for a =0.0 



Figare 5. -Determination of 
proper values 
-J for type of support corre- 
2U spending to 3^"- Eulerian 
2 case as abscissa of curves 
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Pigs. 4,5,6 
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1. Upper oscillation 
Funderaental " 



0, Bulerian load 



a 



1.00 1.50 



Figure 4. -XI for the fundemental 

oscillation and first 
upper oscillation of a bar con- 
formalDly to Suler case, plotted 
against a . 
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m, X = Z 

n, X =s/4 2 

o, X =1/2 Z 

p, Eulerian load 

(point of discontinuity) 

Figure 5. -Deviations of "bar 

form from straight 
line according to (46) for 
x = I ,x = 3/41, and 
X = 1/2 I vertius. a. 
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Pigs. 7,8,9 



e, Eulerian load 
(point of dis- 
continaity) 



Figure T.-Heduction factor - versus a. 
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ra, X = I 

n, X =?/4 1 

o, X =1/2 I 

p, Enlerian load 
(point of dis- 
continuity) 



Figure 8. -Deviations of tar from straight line (according to (49) for 
X = I , X =3/4 1 , and x =l/ S I versus a. 



1.50( 




a 



Figure 9. -Equation (5?). 
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Pigs. 10,11,12,13 




1.00 



Figure 10. -Ratio of fundemental 

frequency for a^iPj; and 



versus ag. 
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Figure 11.- -nl/l-a versus a 



d^n 
m, max £ 

dx2 

n, min " " 




-Eulerian load 
K load 



1.00 



Figure 12.-Uaximun and minimum 

(absolute) of versus a, 

d?2~ 




Figure 13. -Elastic line 

of a "bar 
supported according to 
S^^ Eulerian load case 
and under force P with 
a given eccentricity. 
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Pigs. 14,15 




Pigure 14. -Ratio c of absolute asstimed maxim-um moments in static and 

dynamic case for support according to 3^^ Eulerian case 
versus "b proportional to shock period. (Parameter ajj denotes ratio of 
shock force to Sulerian "buckling force) . 




= 1.0 
= 0.75 
= 0.50 
= 0.35 
= 0.10 



Pigure 15.-Replotting of fig. 14 (c versus ajj, parameter h)^ 
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Figs. 16,17 
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Figure 16 . -Maocimum possible c(for constant shock force and support ac- 
cording to Ealer case)and c values(shock period)versus ajj. 
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Figure 17. -Shock period ("b) for different (aj?) at which static and dy- 
namic stress are equal. 
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Pigs. 18,19 
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Figure 18. -Ratio c for ag = 0.77 and a greater range of shock period (h) . 




m, 1.09(1. 157) 
11,1.05(1.085) 
0,1.01(1.015) 
P.O. 99(0. 984) 
q, 0.95(0. 920) 
r,0. 91(0. 864) 
8,0.84(0.770) 
t, 0.75(0. 667) 
u, 0.50(0.414) 
v,a;=0.25(a5;= 
0.198) 



Figure 19.- 



versTis "b. 
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Pigs. 20,21,22 
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Zigiire 20.-— 
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max 



versus a^, parameter "b. 
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Figure 21. -Elastic line of "bar hinged at cither end and stressed in 
conipression F and single transverse loeul 1. 
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Figure 22. -Elastic line of "bar clamped at one and hinged at the other 
end under load P and a single transverse load 1. 



